Vector theory of gravity in Minkowski space-time: flat Universe without black holes 
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We propose a new classical theory of gravity which is based on the principle of equivalence and 
assumption that gravity, similarly to electrodynamics, is described by a vector field in Minkowski 
space-time. We show that such assumptions yield a unique theory of gravity; it passes all available 
tests and free of singularities such as black holes. In the present theory, gravity is described by four 
equations which have, e.g., exact analytical solution for arbitrary static field. For cosmology our 
equations give essentially the same evolution of the Universe as general relativity. Predictions of our 
theory can be tested within next few years making more accurate measurement of the time delay 
of radar signal traveling near the Sun or by resolving the supermassive object at the center of our 
Galaxy with VLB A. If general relativity is correct we must see a steady shadow from a black hole at 
the Galactic center. If the present theory is right then likely the shadow will appear and disappear 
periodically with a period of about 20 min as we predicted in JCAP 10 (2007) 018. Observation of 
such oscillations will also provide evidence for dark matter axion with mass in meV range. 



I. INTRODUCTION 

In 1915 Albert Einstein [l| completed the general the- 
ory of relativity which then became an accepted theory 
of gravity. In general relativity the space-time geometry 
Qik (metric tensor) is the gravitational field described by 
the action 




(1) 



where G is the gravitational constant and c is the speed 
of light. The second term in Eq. ^ describes interaction 
between gravitational field and matter with the rest mass 
density p{t,r). Variation of (P) with respect to git yields 
Einstein equations 



Rij 



(2) 



where Rik is the Ricci tensor and Tik is the energy- 
momentum tensor of matter. 

Einstein equations ^ are a consequence of the postu- 
late that space-time geometry gik is gravitational field, 
and thus we live in curved space-time. From our point of 
view this postulate is counterintuitive because the well- 
tested Standard Model of particle physics (which includes 
the electroweak theory and quantum chromodynamics) is 
a field theory in Minkowski space-time. We believe it is 
unlikely that the fourth fundamental interaction, gravity, 
should be described by a theory which such dramatically 
different from the Standard model as general relativity. 
One should mention that so far general relativity was 
tested only at weak gravitational field [2] and thus it is 
not a theory fully confirmed experimentally. 

General relativity also predicts existence of singulari- 
ties such as black holes when a massive star collapses into 
a point with zero volume and infinite matter density. One 
can argue that general relativity becomes invalid in the 



vicinity of singularities and a quantum theory of gravity 
will remove them. In contrast, the present theory is free 
of such singularities at the classical level. 

Motivated by the well-tested Standard model, here we 
propose a new classical theory of gravity which is a field 
theory in Minkowski space-time and based on the princi- 
ple of equivalence. We postulate that space-time we live 
is flat Minkowski space-time and matter does not affect 
space-time geometry. Also we assume, similarly to elec- 
trodynamics, that gravitational field is a 4— vector Ak 
which lives in Minkowski space-time. 

Next we derive action for Ak- The principle of equiv- 
alence states that motion of test particles in the fixed 
gravitational field Ai^ in Minkowski space-time is equiv- 
alent to motion in curved space-time with a metric fik 
which is determined by the field Ak- Thus, the equivalent 
metric fik must be a functional of the vector field Ak ■ In 
Appendix A we show that the principle of equivalence 
gives the following unique answer for fik 



a2 2AiAk . , / .2\ 
fik = ri^ke - 



A^ 



■sinh(^^), 



(3) 



where r^ik 
sor and 



=diag(l, —1, —1, —1) is Minkowski metric ten- 



A^ = il'^'AAk 



(4) 



General covariance (invariance under general coordi- 
nate transformations) is a mathematical device used to 
implement the principle of equivalence 0]. Thus, action 
for gravitational field written in terms of the equivalent 
metric fik must have general covariant form given by 
-^GRi where /qr is the action of general relativity ([T]) in 
which gik is replaced by fik- As a result, we obtain the 
following expression for the action 



I = -- 



167rG 



-ff'^Rrk - 



dx^ dx'^ 4 



(5) 
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where / = det{ fik) = —e^^ and is the tensor in- 
verse to Jik [hm]^™ = <5f )• The first term in Eq. ^ is 
Einstein-Hilbert action in which the Ricci tensor is 
formed from the field /^fc. 

Action ([ni) is a functional of the gravitational field ■ 
Variation of ^ with respect to Ak yields four equations 
for gravitational field. In general relativity all compo- 
nents of are treated as independent under variation 
of the action ([5]). In the present theory this is not the 
case due to constraint imposed by Eq. ([3]). 

Action ([5]) is written in Minkowski metric in covariant 
form, has no free parameters and serves as a foundation of 
the present theory of gravity. Our derivation of the action 
([5]) is unique and, hence, the vector theory of gravity 
which obeys the principle of equivalence is also unique.^ 

In Section IVII we show that current theory passes all 
available tests. At strong field our theory substantially 
deviates from general relativity and yields no black holes. 
For cosmology the present theory gives essentially the 
same evolution of the Universe as general relativity. 

A remarkable feature of our theory is that equations for 
gravitational field can be solved analytically for arbitrary 
static mass distribution (see Sec. lIVp . If point masses 
are located at ri, r2, . . . then exact solution is 



fik 



where 
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(7) 



and mfe (fc — 1, . . . , N) are constants determined by the 
value of masses. 

Solution ^ is free of black holes. For a star of mass 
M and radius R Eq. ^ reduces to 0(r) = —GM/c^r 
(r > R) and using Eq. ([46]) for energy conservation we 
obtain that escape velocity for a particle from the stellar 
surface is 



V = CsVl-e2'?^(«), 



(8) 



where Cs — ce^"^*^^^ is the speed of light at the stellar sur- 
face (see Eq. (|5l| ). Eq. ([5]) shows that escape velocity 



^ One should mention that, similarly to the vector potential in 
electrodynamics, ylj. is not observable directly. Hence, Aj, can 
be any vector field which gives physically reasonable equivalent 
metric /^j,. Because /^^ depends on Af^ quadratically the equiv- 
alent metric ((Sjl remains real also for pure imaginary A^- If A^ 
is real (pure imaginary) then always /oo < 1 {/oo > !)• We will 
allow both real and pure imaginary Af^ . If gravitational field re- 
alizable in nature corresponds only to real A^. , then solutions of 
our equations for particular physical problems will automatically 
yield real A^. This is, e.g., the case for static gravitational field 
produced by positive masses (see Eqs. (|6} and (0 which yield 
/oo < 1). 



is always smaller then Cg (Cg < c). In addition, solution 
© predicts that starts do not collapse into a point sin- 
gularity but rather form stable compact objects with no 
event horizon and finite gravitational redshift 

In recent years, the evidence for the existence of ultra- 
compact supermassive objects at centers of galaxies has 
become very strong. It is important to note that present 
solution ([6]) not only argues that such objects are not 
black holes, but also can explain quantitatively their ob- 
served properties (see Sec. IVIIII and Ref. (5j). 



II. EQUATIONS FOR GRAVITATIONAL FIELD 
IN MINKOWSKI SPACE-TIME 



Here we obtain equations for gravitational field from 
the action ([5]). In the rest part of the paper raising and 
lowering of indexes is carried out using Minkowski tensor 
ry**^ =diag(l, — 1, — 1, — 1). It is convenient to introduce 
new independent functions, a scalar 



' 2 



(9) 



and replace A^ by a unit vector which we also will denote 
as Au 



Ak 



Ak 



(10) 



The new vector Ak obeys the normalization constraint 

AkAf" = 1. (11) 

Because expression under the square root in Eq. (fTU]) 
can be both positive or negative the unit vector Ak is a 
real or pure imaginary vector. Observable quantities, of 
course, are always real. 

We treat spatial components Aa (a = 1, 2, 3) as inde- 
pendent, while 



Ao - ^l + Aj + Al + Al 



(12) 



In terms of new independent functions, 4> and the unit 
vector Ak, the equivalent metric ([3]) has the form 



(13) 



f^k = n.ke-^^ + 2AAk sinh(20), ^ - g-^^ 
and the inverse tensor is 

^ ^^feg20 _ 2A'A'' sinh(20). (14) 
One can find variation of the action ([S]) using formulas 



-fimfklSf 



ml 
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and expressing variation of 5/**^ in terms of d(j) and 5 A". 
Variation of ([5]) with respect to </> and A"' yields the fol- 
lowing four equations in Minkowski space-time 



RikA = 



SttG 



Tik — -^fikT ] A , 



(15) 



where the Ricci tensor i?^^ and Christoffel symbols F'^, 
are formed from the equivalent metric fik 



n _ ik 
^ik — ^ I 

ox'- 

1 
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V dx 



k 



dx'- 



dx^ 



(16) 



(17) 



Tik is the energy-momentum tensor of matter 

T.,k^^e^^hafkbu^u\ T = r''T,k^^e't (18) 
7 7 

p is the rest mass density, is the particle 4— velocity 



(19) 



V" = dx^/dt {a — 1,2,3) is three dimensional velocity 
of particle {V^ = 



Vi), and 



^J (c2 - \/2) g-20 _^ 2 (Aoc + AaV^f sinh(20) 

(20) 

is the 7— factor in the presence of gravitational field. 

Covariant Eqs. (fTSj) for the scalar cj) and the unit vector 
Ak are the main equations of the present theory of grav- 
ity. They are written in Minkowski metric which means 
that raising and lowering of indexes is carried out us- 
ing Minkowski tensor. In our theory motion of particles 
in gravitational field is described by Eq. (jBSP identical 
to those in general relativity in metric fik- We obtain 
equations of particle motion in Appendix [B) 

One should note that because (j) and Ak are not ob- 
servable directly they may not be smooth functions and 
components of Ak can have jumps at points where = 0. 
However, equivalent metric fik constructed from and 
Ak must be continuous and smooth. 

Present Eqs. have exact analytical solutions in 

the class of functions for which only one component of 
Ak (e.g. component m) is nonzero and <j) is independent 
of x™. In particular, solution Ak = [1,0,0,0], (j) — 0(r) 
describes static field, while solution with Aa — i corre- 
sponds to a gravitational wave with polarization along 
the a— axis. We discuss such solutions next. 



III. GRAVITATIONAL WAVES OF ARBITRARY 
AMPLITUDE 

In general relativity there exist transverse gravitational 
waves of two polarizations 6]. Here we show that the 



same situation takes place in the present theory of grav- 
ity. Let us consider a plane gravitational wave in free 
space. Then gravitational field obeys equations 



R,.kA'' = 0. 



(21) 



Exact solution of Eqs. pTjl for arbitrary strength of grav- 
itational field is given by 



^fc = [0,0,i,0], ^ = ^{t,x,z) 



(22) 



for a wave of y polarization propagating in the x~z plane 
and 



^fc = [0,0,0, i], (f> = (f>it,x,y) 



(23) 



for a wave of z polarization propagating in the x — y 
plane. For both polarizations, is a function satisfying 
the linear wave equation 



□0 = 0, where □ = A - 



(24) 



that is wave propagates with the speed of light c. 

For example, for y polarization (|22p the equivalent 
metric (1131) has the form 





e-20 





































V 














fik — 



and the corresponding Ricci tensor is 



(25) 



R,k = 



/-I 

10 

-e'*'^ 

V 1 



- 2 



(26) 



Obviously Eqs. pT|) are satisfied provided 0(i, x, z) obeys 
the wave equation (j24p . Combining y and z polarizations 
one can find a general solution for arbitrary plane wave 
propagating along the cc— axis: 

Afc = «[0, 0, sina, cos a], (j) = 4>{t, x), a — a{t,x). 

(27) 

The corresponding equivalent metric is given by 



M = 6-2-^7?,^ - sinh(20) 



/O 



2sin^a sin(2a) 

\ sin(2a) 2 cos^ a J 



(28) 

For the gravitational field Eqs. lead to the 

following equations for and a 



n(j) + sinh(4(?!') 



1 



da 
'dt 



da 
dx 



(29) 



,' da d(j) 1 dad(j)\ 
smh(20)n« + 4cosh(20)(-^---^Uo 



(30) 
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which have exact analytical solution in the form 
(j) = 4>{x — ct), a = a{x — ct). 



(31) 



Eqs. and ([?T|) give a general solution for gravita- 

tional wave of arbitrary amplitude propagating along the 
positive a;— direction. 

One should note that in the weak field limit solutions 
obtained above reduce to those of general relativity. In- 
deed, in the weak field limit one can omit nonlinear terms 
in Eq. ()26p . The linearized Rik satisfies Einstein equa- 
tions in free space Rik = and, thus, the present solu- 
tions are also solutions in general relativity. However for 
strong field this is not the case. 

In Appendix C we show that in our theory radiation 
of weak gravitational waves is given by the same formula 
as in general relativity. However, the present vector the- 
ory of gravity is not equivalent to general relativity in 
the weak field limit and there are situations when two 
theories give different answers even for weak field. As an 
example, let us consider solution of Eqs. ((29)) and ([30]) 
in the form 



a = ujt, (j) ~ (j){x), 



(32) 



where u; is a constant. Then Eqs. ([29)1 and ((30)) yield the 
following equation for 0(a;) 



Indeed, for metric ((28)) with a = ckt and 20 = 
C sin(2fca; -I- 6) the Ricci tensor in the first order in (p 
reads 



10 
0-100 





(39) 



Because i?oo a-nd Rn do not vanish our solution ()36p dis- 
obeys Einstein equations ((38)) even for weak field. How- 
ever, present equations (|2ip 

i?22 sin a + i?23 cos a = 0, 



i?33 cos a + i?23 sin a = 

are automatically satisfied because they do not involve 
Rqo and 

Such an example demonstrates that our theory of grav- 
ity is not equivalent to general relativity in the weak field 
limit (see also Appendix D). This is expected because a 
vector theory can not be equivalent to a tensor theory. 
However, in the weak field regimes for which general rela- 
tivity was tested (such as radiation of weak gravitational 
waves and weak static field) both theories give the same 
answer. 



o2 1 

— I + A:^ sinh(40) = 0, k^- (33) 
ox c 

which has exact analytical solution in terms of the Jacobi 
elliptic function sn(M, m) 



3inh(2(/)) = Csn {2kx + 6, -C^) 



(34) 



where C and 6 are arbitrary constants. 

In the weak field limit {\(f>\ <C 1) Eq. ((55)) reduces to a 
simple harmonic oscillator equation 



dx^ 



4/c^ 







and has a solution 



(35) 



(36) 



Solution given by Eqs. (|32p. ((34)) and ((36)) describes a 
standing helical gravitational wave for which fik = tjik 
at nodal points and the unit vector 



2(f> = Csm{2kx + i 



Ak — i[0, 0, sin(aji), cos(wt)] 



(37) 



rotates in the y — z plane with angular frequency u>. 

Even in the weak field limit solution ((55)) is not a so- 
lution of Einstein equations 



R^k = 0. 



(38) 



IV. STATIC GRAVITATIONAL FIELD 

Here we consider gravitational field produced by rest 
matter with density p{r). In this case the energy- 
momentum tensor ((18p has only one nonzero component 
Too = pc^e^'^ and Eqs. ((T5)) have the following exact 
analytical solution 



Ak = [1,0,0,0], 
and 0(r) obeys the equation 

47rG ^ 



A</> 



(40) 



(41) 



For solution ()40p the equivalent metric is given by 
( e2'^('-) 

fik ' 



V 



_p-20(r) 






-e-20(r) 





,-20(r) 



(42) 



and the corresponding Ricci tensor is 
Rik = 



/ e*-^ 

10 

10 

V 1 



dx^ dx^ 



(43) 



Because Rao = Eqs. ((T5)) are automatically satisfied 
for i = 1,2,3. 
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In Newtonian limit Eq. (|4ip reduces to A(/) = AirGp/c^ 
and, thus, c^0(r) has a meaning of gravitational poten- 
tial. 

Solution (I42p is free of black holes for any mass distri- 
bution and field strength. For a point mass M located 
at r = Eq. jH]) leads to c^A^ = AtiGM5{y) and has a 
solution = —GM I (?r . For N point masses at ri, . . . , 
Eq. dUD yields 

Ac/) = 47r [mi^(ri) + . . . + TOjv5(rjv)] , (44) 

where mi, . . . , ttiat are positive constants. Solution of 
Eq. dSl) is 



r = 



mi 



mAT 



(45) 



Next we consider motion of a particle with rest mass 
m in static gravitational field Equation of particle 

motion in general case is obtained in Appendix B. Eq. 
(lB5t for field (HI) reduces to 



d(e^'^7) 
'dt 



d (7e-2'^V) 



dt 



-7c 



= 20 , 



-20 



(46) 



(47) 



where V0 d(t>/dr^ r = x", V = dr/dt is the particle 
velocity and 



7 = 



(48) 



One can also find equation of particle motion (j47ll di- 
rectly from Lagrange's equation = §7-, where the 
Lagrangian (jB7|) for static gravitational field reads 



=20 



-20 



(49) 



Eq. ((46|) follows from Eq. ([47)1 if multiply both sides of 
Eq. (IT71) by 7e~^'^V and make simple algebraic transfor- 
mations. 

Lagrangian (|49p gives the following expression for the 



particle generalized momentum p = ^ 

p = 7e^^'*mV, 
and particle Hamiltonian H — L 



(50) 



H = e'^'t'jmC^ = ^^2c4g20 _^p2p2g40^ (51) 

Thus, Eq. (|46p is the equation of energy conservation 
E =const, where 



E ~ e^'^7mc^ 



e^mc 



(52) 



is the particle energy and Eq. (|47|) is the equation for 
momentum. 

For a massless particle one should use Eq. (|B9[) which 
for a static field reads 



(53) 



Eq. ([53]) describes propagation of a massless particle with 
speed 



o20 



(54) 



One can see that speed of light depends on the gravita- 
tional field (p and t; < c if is given by Eq. (|45p with 
positive masses. By proper rescaling of coordinates in 
Eq. ([55)1 one can remove the factor e"'*'^ at any given 
point. Let us fix = at infinite distance from masses. 
If an observer at infinity sends a light signal towards the 
Sun then near the solar surface < and light will 
propagate with a smaller speed. This is the explanation 
of Shapiro time delay in the present theory of gravity. In 
our theory the space-time geometry is fixed everywhere 
and given by Minkowski metric. Light signal traveling 
the same distance arrives with a delay if the light trajec- 
tory passes near the Sun. The delay occurs because the 
speed of light is smaller near the solar surface. 

Since Eq. does not contain t explicitly the pho- 

ton frequency luq (measured in time t) remains the same 
during light propagation. However, physical processes 
occur with different rates at different (f). Gravitational 
field (j42p can be removed at a given point by rescaling 
time in the factor V/oo = e"^ (t = r/e"*) and spatial co- 
ordinates by \/—faa = e^"^. In such rescaled coordinates 
identical atoms emit light with equal frequencies lo oc 
dx/dr = e'~'^dx/dt. Thus we obtain 



(55) 



where luq is the photon frequency measured in time t. 

Eq. ([55]) shows that if light emitted by an atom prop- 
agates into a region with larger gravitational potential 
then the detected light frequency is smaller then those 
an identical atom emits at the detection point. This 
phenomenon is known as gravitational redshift of light. 
Eq. ((55|) also shows that in our theory there are no black 
holes. Indeed for the gravitational field created by a point 
mass M: (f> = —GAI/c^r. Therefore if a photon is emit- 
ted at a distance r from the mass M with frequency uj 
then an observer at infinity will detect the photon with 
the energy 



hujQ ~ fkue 



-GM/c^r 



(56) 



According to Eq. ([5S|) no matter how close the photon 
is emitted to the mass M the photon's energy at infinity 
never becomes zero. This means that photon can escape 
from the mass M from any distance. Such a conclusion 
is dramatically different from prediction of general rela- 
tivity. In Einstein's theory photons become trapped by 
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the mass M if they are emitted from a distance smaller 
then the event horizon (that is point mass M behaves as 
a black hole) . 

In Section IVll we show that our theory passes all tests 
of general relativity. In particular, Eqs. ([46]) and ([47]) 
explain correctly the precession of the perihelion of Mer- 
cury and Eq. the gravitational redshift of light. 



V. STATIONARY GRAVITATIONAL FIELD 

Next we consider gravitational field produced by sta- 
tionary mass currents, so that the energy-momentum 
tensor of matter is independent of time. We assume 
that matter velocity V is much smaller then the speed 
of light. However, the scalar potential is not neces- 
sarily small. In this case one can look for solution for 
the equivalent metric (|13p in the form (in the Cartesian 
coordinate system) 



where 



,20 



fik 



\~A. 



36 



20 



Aie2 
— e~' 





-20 






_p-20 







(57) 



where the three dimensional vector A = {A\^ Ai^ Aj^ is 
small, A <C 1, and = 0(r) is arbitrary. For metric ((57)) 
the Ricci tensor upto terms linear in A reads 



Roa 



curl(e'''^curlA) 



Ra 



and Eqs. (fTS]) reduce to 



'a/3 



A0 



eV, V(/)-A = 0, 



curl 



(e'''^curlA) 



"^-^curlAl = ^efpV, 



(58) 



(59) 



here V —{V^, V'^, V^) and p are the matter velocity and 
density respectively. Eqs. (|57)) . ([58]) and ((59|) are valid 
in the first order in V/c and for A = reduce to those 
for a static field. 

For a point mass M rotating around the z— axis with 
an angular momentum L Eqs. ([57)) . ([55]) and ([55]) give 
the following expression for fik in spherical coordinates 
{ct, r, 6, ip) 



fik — 



( e"^ 



y Ar sin i 









Ar sin 6 


e-20 











_j.2g-20 











-r^ sin^(6')e-20 



(60) 



GM 



A 



3c Lr sin 9 



-20 



.20 



' 16M3G2 L V ^ 

(61) 

In the weak field limit GM < c^r Eq. ([6l|) yields 
yl = 2GLsm{9)/c'^r'^ which coincides with the answer 
obtained in general relativity. 

One can find exact analytical solutions for gravita- 
tional field of any strength by making Lorentz transfor- 
mation of Eq. (|42|) . For example, gravitational field pro- 
duced by a relativistic point particle moving with con- 
stant velocity V along the a;— axis is given by 

Ak^ ^=L==[c, -V,Q,0], 



1)] 



and 



Vc2 - y2 ' 

GM 



(62) 



(63) 



where M is the particle rest mass. 

VI. TESTS OF THE THEORY OF GRAVITY 

According to the present theory we live in Minkowski 
space-time in which there is vector gravitational field. 
However descriptions of particle motion in terms of the 
vector field or space-time geometry are equivalent. This 
is Einstein equivalence principle. One can consider in- 
teraction of particles with the field in two different, but 
equivalent ways. In the first approach the space-time 
is Minkowski space-time in which there is vector gravi- 
tational field and motion of a test particle is described 
by Eq. (jB5p . In the second treatment we describe in- 
teraction of the particle with the gravitational field in a 
geometrical way. Namely we assume there is no vector 
field, but instead the space-time is curved with metric 



„ f „ik rik 

gik — jik, g —J , 

and particles move along geodesic lines 
ds 



(64) 



(65) 



\/ gikdx^dx'^ and 

— dx^ /ds is the particle 4— velocity. Please note that 
B5[) is written in Minkowski metric, while Eq. 



where F^^, are Christoffel symbols, ds 
u 

Eq. ^ ^ 
is written in metric (|64p . However, mathematically equa- 
tions (jB5P and (p5|) are identical and thus give the same 
answer. 

It is convenient to compare the present theory of grav- 
ity with observations based on the equivalent metric. Ac- 
cording to Eq. (132) the equivalent metric for static grav- 
itational field is 



gik — 



20 














_e-20 














_e-20 














„e-20 



(66) 
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where for a point mass M: 4> = —Mjr (here we put 
G = c = 1 y. Metric (pS]) is known as Yilmaz exponential 
metric 0, Q • For a point mass M Einstein equations ([2|) 
yield Schwarzschild metric which in isotropic Cartesian 
coordinates has the form [o| 



where 



h{r){dx'y - 9{r)[{dx'y + {dx^ + [dx^], (67) 



h{r) 



1 - M/2r 
1 + M/2r 



9{r) 



M 
27 



(68) 



For small M/r both the Schwarzschild 
(I66p metrics yield the same expansion 



and Yilmaz 



2M 2A'P , , 2M , , 

/i(r) = l + +..., g(r) = l + + ...(69) 



which is known as Post-Newtonian approximation. The 
four classic tests of general relativity, namely the grav- 
itational redshift of light, the deflection of light by the 
Sun, the precession of the perihelion of Mercury and time 
delay of a radar signal traveling near the Sun (Shapiro 
delay) , have examined the metric in the Post-Newtonian 
approximation (j69p 0. Because the equivalent metric 
([66|l obtained in the present theory has correct Post- 
Newtonian limit ([M]) our theory of gravity also passes 
the four classic tests. For static field the present theory 
of gravity gives answer different from general relativity in 
the next correction beyond the Post-Newtonian approxi- 
mation. So far, however, gravity have not been tested in 
this region. 

In Appendix C we show that in our theory radia- 
tion of weak gravitational waves is described by the 
same formula as in general relativity. Such radiation 
was indirectly detected as energy loss by binary pulsars 
and served as a quantitative test of Einstein equations 
for weak time-dependent field. The present theory also 
passes this test. 



VII. COSMOLOGY 

In this section we apply our theory to evolution of the 
Universe. We assume that one can omit kinetic energy of 
matter compared to its rest energy. We also assume that 
matter is uniformly distributed in space with density p, 
where p is independent of time. According to the present 
theory we live in Minkowski space-time and galaxies lo- 
cated in different parts of the Universe do not move rela- 
tive to each other (apart from local random motion which 
we omit in this section). That's why p is constant. How- 
ever gravitational field evolves with time. Light emitted 
by a distant source reaches an observer on Earth with a 
delay. As a result, at the moment of light detection the 
gravitational field is different from its value when light 
was emitted. This is the origin of cosmological redshift. 



For spatially isotropic Universe gravitational field must 
have the form 



Afc = [1,0,0,0] 
and, thus, equivalent metric (jl3p is given by 

/ ^ \ 

, _ -a^ 

-a2 

\ y 



(70) 



(71) 



where a = e~"^. As in general relativity, we add the 
cosmological constant term to energy-momentum tensor, 
Tih Tik + Afik, where A is a constant. Then for grav- 
itational field Eqs. ^ read 



RiO 



SttG 



T, 



iO 



■^fioT — A/io 



(72) 



We also assume that the Universe is homogeneous on 
large spatial scales and, therefore, a in Eq. ((7T|) must 
be independent of r. Then for the equivalent metric ([7T|) 
the Ricci tensor is 



3 d 

Roo = ttt: (aa) , RaO = 0, 

ot 



and 



R. 



■af3 



6c.f3-g-^ [a a) 



(73) 



(74) 



The energy-momentum tensor of matter (fT8|) has only 
one nonzero component Tqq — pc^/a^, T = pc^/a^, and 
therefore Eqs. ((72)) reduce to 



3- (aa) 



AnG ( pc 



2A 



(75) 



Multiplying both sides of Eq. ([TS]) by ad and integrating 
over time we obtain 



3a^ 



SttG f pc 



A- 



(76) 



where G is an integration constant. For A = sign of G 
determines whether a Universe is open or closed. In our 
theory the constant G produces the same effect on a{€) 
as the curvature of space in general relativity. 

In the metric ([71]) Einstein equations ([2]) with the cos- 
mological constant term give our Eq. (|76|) with G = 0. 
Indeed, present Eqs. ([7^ coincide with Einstein equa- 
tions ([2]) with index fc = 0. In addition, Einstein equa- 
tions with ik = aa yield 



9 . o . . SttG f pc^ , , 

— (a'^a) = [ - h Aa^ 

dt^ ' \2a 



(77) 



Multiplying both sides of Eq. (|77)) by a^d and integrating 
over time we obtain 



3a^ = 



SttG 



pc^ , Ci 
^ + A+ ^ 



(78) 
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where Ci is an integration constant. Eqs. ([76)1 and (l78l) 
are compatible ii Ci — C — 0, and thus our Eq. ([75]) 
with C = is also solution in general relativity. Such a 
solution agrees with available cosmological data. 
By changing time coordinate into 

Eq. (Uni) yields 

which describes evolution of the Universe in the metric 
ds"^ = dr'^ — {dx^ + dy'^ + dz^) . Such form of metric is 
commonly used in cosmology for spatially flat Universe. 

VIII. GALACTIC CENTERS AND DARK 
MATTER PROBLEM 

In the present theory static gravitational field is de- 
scribed by the equivalent exponential metric (|66p. Metric 
was also obtained in Refs. 0,[l[i,[ll- Exponential 
metric (1661) predicts no black holes, but rather compact 
objects with no event horizon and very large, but finite, 
gravitational redshift. 

In recent years, the evidence for the existence of an 
ultra-compact concentration of dark mass at centers of 
galaxies has become very strong. However, a proof that 
such objects are black holes rather then compact objects 
without event horizon is lacking. If the present theory 
of gravity is correct then the compact supermassive ob- 
jects at galactic centers can not be composed of baryonic 
matter. Indeed, mass of a compact (neutron star like) 
baryonic object in the exponential metric (j66p can not 
exceed about I2M0 but the objects at galactic cen- 
ters possess masses upto a few IO^Mq. Hence, those 
objects must be made of dark matter of non baryonic 
origin. This fact gives us an opportunity to determine 
composition of dark matter based on observations of su- 
permassive objects at galactic centers. 

In the previous paper Q we showed that properties 
of compact objects at galactic centers can be explained 
quantitatively assuming they are made of dark matter 
axions and the axion mass is about 0.6 meV. Analysis 
of Ref. is based on the assumption that static gravi- 
tational field is described by the exponential metric ([SS]) 
rather then by general relativity. A full time-dependent 
theory of gravity was unnecessary for calculations made 
in Ref. . The present paper provides such a theory and 
justifies our previous choice of the exponential metric. 

Axions are one of the leading particle candidates for 
the cold dark matter in the Universe Interaction of 
axions with QCD instantons generates the axion mass ni 
and periodic interaction potential 

V{^p)=m^F'^[l-cos{Lp/F)], (81) 



where is a real scalar axion field and F is the Peccei- 
Quinn symmetry breaking scale. The interaction poten- 
tial ([5T|l has degenerate minima V = bX ip = 27TnF, 
where n is an integer number. As a consequence, ax- 
ions can form bubbles. Bubble mass is concentrated in 
a thin surface (interface between two degenerate vacuum 
states). In the exponential metric the potential energy 
of a spherical bubble with radius R is given by 5] 

U{R) = inaR^ exp , (82) 

where a is the surface energy density and M is the fixed 
bubble mass. U{R) has a shape of a well. At i? 3> M 
one can omit gravity and U (i?) ~ AnaR"^ is just a surface 
energy (tension) which tends to contract the bubble. At 
R M gravity effectively produces large repulsive po- 
tential which forces the bubble to expand. As a result, 
the bubble radius R{t) oscillates between two turning 
points. 

In Ref. 0, based on quantitative analysis of avail- 
able data, we argued that such oscillating axion bubbles, 
rather then supermassive black holes, could be present at 
galactic centers. Recent observations of near-infrared and 
X-ray flares from Sagittarius A* , which is believed to be a 
3.6 X IO^Mq black hole at the Galactic center, show that 
the source exhibits about 20-minute periodic variability 
[13I . [13, [1^ . An oscillating axion bubble can explain such 
variability. Known value of the bubble mass at the center 
of our Galaxy and its oscillation period yields the axion 
mass of about 0.6 meV. Size of the axion bubble at the 
center of the Milky Way oscillates between i?niin ~ l-R© 
and i?,nax « lAU « 21Oi?0. 

Further, as shown in Ref. the axion bubbles with 
no free parameters (if we fix m = 0.6 meV based on Sagit- 
tarius A* flare variability) quantitatively explain the up- 
per limit (a few IO^Mq) on the supermassive "black hole" 
mass found in recent analysis of the measured mass dis- 
tribution [1^. Also, with no free parameters the bubble 
scenario explains observed lack of supermassive "black 
holes" with mass M < lO^M© [l^. For such low-mass 
bubbles the decay time t oc M^/^ becomes much shorter 
then the age of the Universe and, as a result, such objects 
are very rare. 

Observation of the Galactic center with very long- 
baseline interferometry within the next few years will 
be capable to test theories of gravitation in the strong 
field limit. Such an observation will allow us to distin- 
guish between the black hole (predicted by general rela- 
tivity) and the oscillating axion bubble scenario. A defin- 
ing characteristic of a black hole is the event horizon. 
To a distant observer, the event horizon casts a rela- 
tively large "shadow" over the background source with 
an apparent diameter of about lOGM/c^ « 8Oi?0 due 
to bending of light. The predicted size of this shadow 
for Sagittarius A* approaches the resolution of current 
radio-interferometers. Hence, there exists a realistic ex- 
pectation of imaging the shadow of a black hole with 
VLB A within the next few years [H, [H, H [SI, [H . If 
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the axion bubble, rather then a black hole, is present at 
the Galactic center, the steady shadow will not be ob- 
served. Instead, the shadow will appear and disappear 
periodically with a period of about 20 min. Discovery 
of periodic appearance of the shadow from the Galactic 
center object will also be a strong evidence for the ax- 
ion nature of dark matter and will lead to an accurate 
prediction of the axion mass. 

One should mention that intrinsic size of Sagittarius 
A* at a wavelength of 1.3 mm was recently determined 
using VLB A [231. The intrinsic diameter of Sagittarius 
A* was found to be < 0.3AU~ 65i?0 which is less than 
the expected apparent size of the event horizon of the 
presumed black hole. Such observation might indicate 
lack of black holes, in agreement with the present theory. 



IX. CONCLUSIONS 



Ref. ^ is based on the present exponential metric ([66]) 
for static gravitational field rather then general relativity. 

Our theory of gravity can be tested in several ways. 
For example, one can examine gravity beyond the Post- 
Newtonian approximation in the solar system by improv- 
ing the accuracy of Shapiro time delay experiment (time 
delay of a radar signal traveling near the Sun). Another 
possibility is to resolve the supermassive object at the 
center of our Galaxy with VLBA. If general relativity is 
correct we must see a steady shadow from a black hole. 
If the present theory is right then likely the shadow will 
appear and disappear periodically with a period of about 
20 min as we predicted in Q. Observation of such oscil- 
lations will also provide evidence for dark matter axion 
with mass in meV range. 

This work was supported by the Ofhce of Naval Re- 
search (Award No. N00014-07-1-1084 and N0001408-1- 
0948). 



Here we propose a new classical theory of gravity which 
is based on the principle of equivalence and assumption 
that, similarly to electrodynamics, gravity is described 
by a vector field in Minkowski space-time. We show that 
present theory is the only possibility that can be obtained 
from these assumptions. Our theory fundamentally dif- 
fers from general relativity which treats space-time geom- 
etry as gravitational field. In the present theory, similarly 
to the Standard Model, matter does not affect geometry 
of flat Minkowski space-time. 

The current vector theory is not equivalent to general 
relativity even in the weak field limit. Nevertheless, our 
theory also passes all available tests and for static field 
in the Post-Newtonian approximation gives the same an- 
swer as general relativity. Beyond the Post-Newtonian 
approximation the present theory gives different result 
and yields no singularities such as black holes. A defin- 
ing characteristic of a black hole is the event horizon. So 
far there were no observations of the event horizon and, 
thus, a proof of black holes existence is lacking. For cos- 
mology our theory predicts essentially the same evolution 
of the Universe as general relativity. 

In the present theory gravitational field is described 
by four equations (fT5|) which can be solved analytically 
for much greater number of problems then ten Einstein 
equations In particular, for arbitrary static mass 

distribution our Eqs. (|15p have exact analytical solution 
(113. 

The present theory, if confirmed, can also lead to a 
break through in the problem of dark matter. Namely, 
the theory predicts that supermassive compact objects at 
galactic centers have non baryonic origin and, thus, yet 
undiscovered dark matter particle is a likely ingredient 
for their composition. As a result, observations of such 
objects can allow us to predict the nature of dark mat- 
ter. In the previous paper Q we showed that properties 
of compact objects at galactic centers can be explained 
quantitatively assuming they are made of dark matter 
axions and the axion mass is about 0.6 meV. Analysis of 



APPENDIX A: DERIVATION OF EQUIVALENT 
METRIC 



Here we obtain expression for the metric fik which is 
equivalent to the vector field Ak using the principle of 
equivalence. For small A^. one can expand fik in pow- 
ers of Ak- In the leading order there are two possible 
terms AiAk and A^Tjik, where r/.k =diag(l, -1, -1, -1) 
is Minkowski metric tensor and A^ = AkA'' is square of 
vector Ak (^' = ri'^'^Ak). To obey available experimen- 
tal tests these terms must enter fik in the combination 
A^rjik — 2AiAk and, therefore, fik has the form 



fik ~ Vtk + A^rjik - 2 AiAk. 
One can write Eq. (|A1|) as 



fik ~ Vik + nim't'k^ 



where 



^A^5T 



2A"'Ak 



(Al) 

(A2) 

(A3) 
21) reads 

(A4) 



In curved space-time with metric gik Eq. 

fik ~ gik + gim4>T^ 

where c^™ is written in metric gik- 

Let us now assume that gravitational field is not small. 
According to the equivalence principle one can describe 
motion of particles as if there is no field, but instead the 
space-time is curved with an equivalent metric gik ■ If we 
change (jjik by a small amount (pik — > 4>ik + o:4)ik, where 
a is a small number, then we can write (0 = (pik) 



ik\ 



^4>) ~ fik{<l>) 



df,k[{l + a)<j>] 



a=0 



da 

From the other hand Eq. (|A4|) yields 

/,fc(0 + a<l>) « fM + f^Mac|>T. 



(A5) 



(A6) 
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Comparing Eqs. (|A5p and (jA6[) we obtain that for any 



dc 



(A7) 



a=0 



Because Eq. (jA7l) must be satisfied for any (/(^^ it is 
essentially a differential equation, rather then a condition 
at one point (a = 0). Solution of Eq. (|A7p yields the 
following unique expression for the metric equivalent to 
the gravitational field 4>ik 



fik = ?7™exp((/)"), 
where the exponent of a tensor is defined as 



(A8) 



exp(C) =ST + cf>T + ^<t>T<t>i + li<t^T<t'p€ + • • • • (A9) 



Substituting Eq. into Eq. 

account that 



exp(-2A'Afc) = 6i 



^2 



and taking into 
(AlO) 



we obtain 

U = v^meMA^^T - 2A"Mfc) = e'^'77,„exp(-2A"Mfc) 



^^.^ 2'^-'^"^'=.sinh(A2). (All) 



In Eq. (jAlip Ak is the gravitational field in metric fik- 
Therefore, Ai = fikA'^. Taking into account Eq. (jAlip 
we find 

A^e-^\,kA\ A^^e^^^'^Ak, (A12) 
and A^ — AkA^ is determined from the equation 

A^e-^" = rfAAk- (A13) 
As a result one can write fik as 



f^k = e^V^k 



A,A 



A^ 



(e^- -1) 



(A14) 



It is convenient to rewrite this equation in a different 
form introducing Ak — AkC"^ Then 



U = e^\.k - 2^ sinh(i2), (A15) 



where 



A^^Tj^'^AAk- 



(A16) 



APPENDIX B: MOTION OF PARTICLES IN 
GRAVITATIONAL FIELD 

Here we find how a test particle with rest mass m 
moves in an external gravitational field fik- Interaction 
of the particle with the field is described by the action 

/fiold-mattcr = -mc J ^ fikdx^dx'', (Bl) 

where the integral is taken along the particle trajectory. 
One can find equation of particle motion varying the ac- 
tion dm]) at fixed f,k 6] 



<5/field-matter = ^ J [u'u'' dsS fik + fikiu" d6x' + dSx" 

(B2) 

where is the particle 4-velocity 

dx' 



ds 



(B3) 



with ds = -y/ fikdx^dx^ . Next we take into account 5 fik = 
{dfik/dx'') Sx'' and integrate the second term by parts 



5I& 



cld-mattcr 



mc 



/ dhk_ _ dfik _ dfu 
\ dx^ dx^ dx^ 



du'' 

^fik—r- 

ds 



dsdx 



Principle of least action (5/fieid-matter = yields the fol- 
lowing equation 



duf" _ 1 
■'"'"S" ~ 2 



dfik _ dfik _ dfii 
dx'- dx^ dx^ 



(B4) 



Multiplying both sides of Eq. (|B4p by tensor inverse to 
fik we find 



ds 2^ 



dfik _ dfik _ dfii 
dx^ dx^ dx'^ 



i k 

u u , 



or 



du» 
ds 



(B5) 



(B6) 



This is equation of motion of a particle in gravitational 
field fik- 

From Eq. (|Bip we obtain the following Lagrangian of 
the particle 



L 



-mc\ fik 



dx^ dx'^ 



dt dt 



(B7) 



Action (jBip and Eq. (jBSp are invalid for massless 
particles. Let us consider a massless scalar field x- In 
the gravitational field fik the action for x reads (in 
Minkowski metric) 



ff 



IV dx* dx 
dxt^ dx^ 



(B8) 
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Variation of Eq. (|B8[) yields the following equation of 
motion for the field x 



For geometrical optics one can write x as x 



(B9) 



IX|e 



where ip (eikonal) has a large value. Substituting this into 
Eq. (|B9p and keeping only the leading term we obtain 
eikonal equation in gravitational field 



/ 



0. 



(BIO) 



APPENDIX C: RADIATION OF 
GRAVITATIONAL WAVES 

Here we consider radiation of weak gravitational waves. 
In general relativity the metric for weak plane gravita- 
tional waves propagating along the x— axis in a properly 
chosen coordinate system reads fd\ 



gik = Vik + 







h22{t,x) h23{t,x) 

h2'i{t,x) -h22{t,x) 



(CI) 



where /i23 and /122 are small perturbations. In linear 
approximation the corresponding Ricci tensor is 



Rik — 7: 



/O 



Uh22 0/123 

Vo 0/^23 -0/122 



(C2) 



To find energy radiated in the a;— direction one should 
find functions /123 and h22 in the wave zone. They are 
determined from Einstein equations ^ 



1, 



'-Uh 



8t:G 



22 



To: 



1 



:322T 



1^, SttG 

- 0/123 — —^T23, 



^□/l22 



SttG 



T: 



33 



^533r 



(C3) 



(C4) 



(C5) 



In the right hand side of Eqs. (|C3|) - (|C5p one can replace 
gik with Minkowski metric rjik- 

In the present theory of gravity for weak gravitational 
waves {\(f>\ ^ 1) propagating along the x— axis 



Ak — i[0, 0, sina, cos a] 
and the equivalent metric (|13p has the form 



fik Vik + 20 



-10 
10 

cos(2a) -sin(2a) 

-sin(2Q;) -cos(2a) 



(C6) 



(C7) 



where (p = (j>{x — ct) and a = a{x — ct). Eqs. (jC6p and 
(|C7|) are equations in Minkowski space time. Let us make 
the following change of coordinates 



t + - 
c Jo 



(f>{q)dq, 



(j}{q)dq. 



(C8) 



(C9) 



Taking into account that under coordinate transforma- 
tion x* x* -I- where are small, tensor fih trans- 
forms as fik fik ~ d£,Jdx'' - d^k/dx' [g] we obtain 
that in new coordinates 



fik = Vik + 20 



/O 



cos(2q;) 

V -sin(2a) 




(CIO) 



which is precisely Eq. (jCip with /i22 = 20cos(2q;) and 
/i23 = — 20sin(2Q;). As a result, for weak gravitational 
waves our equivalent metric fik in the wave zone is equal 
to gik, where gik is solution of Einstein equations (jCSP - 
((C5)) . Indeed, if f^k obeys Eqs. (IC3|-(IC5| then it auto- 
matically satisfies present equations 



RikA — 



SttG 



T, 



ik 



\hkT 1 A" 



(Cll) 



which are linear combinations of Eqs. (|C3|1 - (|C5P with 
coefficients A^ . 

We obtain that general relativity and our theory yield 
the same answer for the equivalent metric of emitted 
weak gravitational waves. Therefore, in both theories the 
energy loss due to emission of weak gravitational waves 
is given by the same formula. 



APPENDIX D: WEAK GRAVITATIONAL FIELD 

Here we obtain equations for weak gravitational field. 
In this case \4>\ 1, however, the unit vector Ak can be 
arbitrary. For weak field the Ricci tensor reduces to 



Ri. 



dx^ dx'' 



ik 



(Dl) 



Taking into account Eq. ([T7| for P'^, and formula F'^ = 
d\ny/^/dx' = -2a0/9a?lg we find that Eqs. ^ in 
the weak field limit yield 



2A'= 



92 



dx^dx^ 



{A.A'cf) 



92 



dx'dx 



J (AkA^cj)) + n{A.,Ak^ 



A,Dcb = ^ ( T,fc 



-Tj.kT ] A\ (D2) 
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where 

Tik is the energy momentum tensor for free particles 

(D3) 



Tik — 



1 



-pC U.iUk, 



T= \ l- 



1/2 



,2 ^C', 



V 

p is the rest mass density and Uk is the particle 4— velocity 



1 



Uk 



1 



y2 



1, 
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(D4) 



The weak field equations (|D2[) are Lorentz covariant and 
nonlinear in A^. Eqs. (jD2l) differ from the weak field 
limit of Einstein equations |6||. This is the case because 
the present vector theory of gravity is not equivalent to 
general relativity (tensor theory) even for weak field. 



If Aa are also small ^ 1 and « 1) one can 

rewrite Eqs. (jD2|) in a form similar to Maxwell's equa- 
tions (for consistency we assume that F <C c and omit 
second order time derivatives) 



where 



divE = At: 



Gp 



curlB 



AirGpY 



laE 

c 9F' 



E = V0, B =curl(A0), 



(D5) 



(D6) 



A^{A^,A^,A^) and Y ^{V^,V^,V^) is the three di- 
mensional velocity of particle. 
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